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ABSTRACT 

We report the final results of our study of the cosmic microwave background (CMB) with the BIMA 
array. Over 1000 hours of observation were dedicated to this project exploring CMB anisotropy on 
scales between 1' and 2' in eighteen 6'. 6 FWHM fields. In the analysis of the CMB power spectrum, 
the visibility data is divided into two bins corresponding to different angular scales. Modeling the 
observed excess power as a flat band of average multipole igff = 5237, we find AT^ = 220li25 AiK2 
at 68% confidence and AT^ > /zK^ with 94.7% confidence. In a second band with average multipole 
of ieff = 8748, we find ATj consistent with zero, and an upper limit 880 /xK^ at 95% confidence. 
An extensive series of tests and supplemental observations with the VLA provide strong evidence 
against systematic errors or radio point sources being the source of the observed excess power. The 
dominant source of anisotropy on these scales is expected to arise from the Sunyaev-Zel'dovich (SZ) 
effect in a population of distant galaxy clusters. If the excess power is due to the SZ effect, we can place 
constraints on the normalization of the matter power spectrum cg = LOS^Q ^g at 68% confidence. The 
distribution of pixel fluxes in the BIMA images are found to be consistent with simulated observations 
of the expected SZ background and rule out instrumental noise or radio sources as the source of the 
observed excess power with similar confidence to the detection of excess power. Follow-up optical 
observations to search for galaxy over-densities anti-correlated with flux in the BIMA images, as 
might be expected from the SZ effect, proved to be inconclusive. 

Subject headings: cosmology: observation - cosmic microwave background - Sunyaev-Zel'dovich effect 



1. INTRODUCTION 

The angular power spectrum of the Cosmic Microwave 
Background (CMB) has been measured with high sig- 
nal to noise on scales from degrees to several arcminutes 
(e.g., Hinshaw et al. 2003, Kuo et al. 2004, Mason et al. 
2003). However, observations of CMB anisotropy on ar- 
cminute scales, where secondary anisotropics such as the 
Sunyaev-Zel'dovich (SZ) effect (Sunyaev & Zel'dovich 
1970) are expected to dominate the primary CMB 
anisotropy (e.g., Gnedin & Jaffe 2001), have not yet 
achieved comparable precision. At these finer angular 
scales, observations distant SZ clusters have the potential 
to be a powerful probe of the growth of structure in the 
Universe (Carlstrom, Holder, & Reese 2002). However, 
the power spectrum of the arcminute fluctuations reveals 
little information about the nature of the sources respon- 
sible for the anisotropy. As suggested by Rubino-Martin 
& Sunyaev (2003), higher order statistics of images can, 
in principle, be used to identify the unique signature of 
the SZ effect. 

Beginning in the summer of 1998, we began a pro- 
gram to search for arcminute-scale CMB anisotropy us- 
ing the Berkeley- Illinois-Mary land Association (BIMA^) 
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interferometer. Initial results are included in earlier pub- 
lications (Holzapfel et al. 2000, Dawson et al. 2001, and 
Dawson et al. 2002 (hereafter D2002)). A detailed de- 
scription of the BIMA analysis, results from other ex- 
periments, and a comparison with theoretical models 
and simulations of structure formation can be found in 
D2002. In this paper, we report the final results from 
the BIMA CMB anisotropy survey. We describe the ob- 
servations of the final fields in summer 2002 with the 
BIMA array and the Very Large Array (VLA^) in 321 
The results of the power spectrum analysis, including a 
discussion of tests for systematic errors in the analysis 
are presented in ^ In ^we examine the BIMA image 
statistics in order to constrain the origin of the observed 
anisotropy. In ^ we describe the results of follow-up op- 
tical observations used in an attempt to identify galaxy 
clusters in the BIMA fields. Finally, in ^ we summarize 
the results and present our conclusions. 

2. OBSERVATIONS 

The BIMA anisotropy survey consists of 18 fields that 
had not been previously observed for SZ galaxy clusters 
or CMB anisotropy at arcminute angular scales. Each 
field was observed with the BIMA array at a frequency of 
28.5 GHz and the VLA at a frequency of 4.8 GHz. Anal- 
ysis of ten fields observed with the BIMA array during 
the summers of 1998, 2000, and 2001 revealed evidence 
for a detection of power in excess of the instrument noise. 
A description of the analysis and results for the first ten 
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fields in the survey can be found in D2002. Eight new 
fields were added to the survey in the summer of 2002. 

2.1. BIMA Observations 

All anisotropy observations were made using the BIMA 
array at Hat Creek. Nine 6.1 meter telescopes of the ar- 
ray were equipped for operation at 28.5 GHz, providing 
a 6.6' FWHM field of view. The first ten fields, BDF4- 
BDF13, were chosen to lie in regions expected to have 
minimal radio point source and dust contamination. The 
eight fields added in 2002, BDF14-BDF21, were chosen 
to lie two minutes east in Right Ascension of existing 
fields in the survey. Observations of these fields added 
another 489 hours of observations to the previously pub- 
lished data, for a total of 1096 hours for the complete 
BIMA survey. Each of these new fields was observed us- 
ing the same phase calibrator as its previously observed 
nearest neighbor and analyzed following the data reduc- 
tion as described in D2002. These fields were chosen to 
check for possible contamination of the signal correlated 
with telescope position. This makes for a total of 18 inde- 
pendent fields in the survey, covering approximately 0.2 
square degrees. The pointing center, dates of observa- 
tion, and observation time for each of the fields are given 
in Table ^ We dedicated 55 — 65 hours of observation 
with the BIMA array to each of the new fields in order 
to achieve a uniform RMS noise level of < 150/zJy/beam 
on short baselines {u-v < 1.1 kA) for the entire sample. 
This noise level corresponds to an RMS of 15.55 fiK for 
a 2' synthesized beam. 

2.2. VLA Observations 

To help constrain the contribution from point sources 
to the anisotropy measurements, we used the VLA at 
a frequency of 4.8 GHz to observe each field in the sur- 
vey. With 1.5 hours per field, these observations yielded 
an RMS flux of ~ 25 ^Jy/beam over a 9' FWHM region 
with the same pointing center as a BIMA field. The posi- 
tions of all point sources detected with significance > 6a 
within 400" of the pointing center have been recorded. 
Measured point sources with fluxes corrected for attenua- 
tion by the primary beam at 4.8 GHz are listed in Tables 
Hand El 

If the spectra of the point sources are nearly flat or 
falling, deep observations with the VLA will identify 
those that lie near the noise level in the 28.5 GHz maps. 
However, it is possible that a radio source with a steeply 
inverted spectrum may lie below the VLA detection 
threshold but would still contribute significantly at 28.5 
GHz. Advection dominated accretion flows are thought 
to be the most common inverted spectrum sources. They 
typically have a slowly rising spectrum, with a spectral 
index of 0.3 to 0.4 (Perna & DiMatteo, 2000) where the 
point source flux S" oc z^" . Such a shallow spectrum would 
only provide a factor of two increase in flux between 4.8 
GHz and 28.5 GHz; any source not seen with the VLA 
would be near the noise level of the BIMA observations. 

To search for point sources with more steeply inverted 
spectral indices, we made VLA observations at 8.0 GHz 
of the five BIMA fields that most strongly indicate an 
excess of anisotropy power. If radio sources are the dom- 
inant contribution to the observed excess power, it is 
these fields that are the most likely to be contaminated. 
Results of these observations are found in Table |21 and 



Table |3| The 8 GHz observations reached a RMS flux 
density of 20 /xJy/beam at the center of the 5' FWHM 
primary beam. Six of the 13 sources identifled in the 
4.8 GHz maps were detected in the 8.0 GHz maps with 
Signal-to-Noise ratio (SNR) > 3. The mean spectral in- 
dex of the detected sources is found to be a = —0.4. 
The fluxes of the other seven sources were poorly con- 
strained because those sources were either too dim or 
were positioned outside of the 8.0 GHz primary beam 
where the instrument is most sensitive. The 8 GHz im- 
ages produced no additional point source detections with 
SNR> 6, providing additional evidence against contam- 
ination by a population of radio sources with inverted 
spectra. 
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TABLE 1: Field Positions and Observation Times 



Fields R. A. (J2000) Peel. (J2000) Observation year(s) Time (Hrs) 



BDF4 




28"" 


04.4" 


+28° 23' 06" 


98 


77.6 


HDF 


12'^ 


36"" 


49.4= 


+62° 12' 58" 


98, 01 


59.9 


BDF6 


18'' 


21"' 


00.0= 


+59° 15' 00" 


98, 00 


81.2 


BDF7 


06'^ 


58'" 


45.0= 


+55° 17' 00" 


98, 00 


68.2 


BDF8 


OO'' 


17'" 


30.0= 


+29° 00' 00" 


00, 01 


53.3 


BDF9 


12'' 


SO'" 


15.0" 


+56° 52' 30" 


00, 01 


53.9 


BDFIO 


18'' 


12"" 


37.2= 


+58° 32' 00" 


00, 01 


53.3 


BDFll 


06'' 


58'" 


00.0= 


+54° 24' 00" 


00, 01 


50.0 


BDF12 


06'' 


57"' 


38.0= 


+55° 32' 00" 


01 


54.8 


BDF13 


22'' 


22"' 


45.0= 


+36° 37' 00" 


01 


54.5 


BDF14 


00'' 


26"' 


04.4= 


+28° 23' 06" 


02 


62.4 


BDF15 


06'' 


56"' 


45.0= 


+55° 17' 00" 


02 


64.2 


BDF16 


12'' 


34"' 


49.4= 


+62° 12' 58" 


02 


64.5 


BDF17 


18'' 


19™ 


00.0= 


+59° 15' 00" 


02 


64.5 


BDF18 


00'' 


15"' 


30.0= 


+29° 00' 00" 


02 


57.8 


BDF19 


06'' 


55"' 


38.0= 


+55° 32' 00" 


02 


59.4 


BDF20 


12'' 


48"' 


15.0= 


+56° 52' 30" 


02 


54.7 


BDF21 


18'' 


10"' 


37.2= 


+58° 32' 00" 


02 


62.0 
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TABLE 2: Point Sources and Fluxes Identified from VLA Observations(1998-2001) 



Field 


A R.A. (") 


A DEC (") 


4.8 GHz (^Jy) 


8 GHz (/iJy) 30 GHz (/xJy) 


BDF4 


-96.8 


255.7 


1230 ±90.6 






BDF4 


72.8 


178.2 


514 ±54.9 




468 ±313 


BDF4 


99.9 


-89.4 


221 ± 42.3 






BDF4 


-94.9 


268.4 


391 ± 89.3 






HDF 


-35.0 


-85.0 


832 ± 56.3 




535 ± 1 85 


HDF 


255.0 


-89.8 


1380 ± 93.0 






HDF 


178.3 


-274.0 


1520 ± 112 






HDF 


222.5 


-86.8 


709 ± 64.6 






HDF 


69.1 


334.1 


1120 ± 107 






HDF 


-21.1 


66.1 


190 ±36.2 






BDF6 


-136.5 


-283.5 


592 ± 74.5 


—51 3 ± 226 




BDF7 


314.6 


47.4 


1554 ± 80.5 






BDF7 


173.8 


97.8 


373 ±40.1 




156 ± 306 


BDF7 


253.8 


- 1.1 


284 ± 49.5 






BDF8 


-145.9 


-266.1 


1381 ±94.1 






BDF8 


27.6 


280.9 


611 ± 72.8 






BDF8 


302.9 


-79.5 


622 ± 86.1 






BDF9 


-221.9 


-123.7 


1500 ± 78.4 






BDF9 


-192.7 


215.8 


1193 ±81.3 






BDF9 


245.2 


-101.0 


1039 ± 67.8 






BDFIO 


-158.5 


-165.6 


1670 ± 63.9 




1610 ± 347 


BDFIO 


-146.1 


-183.9 


320 ± 44.4 






BDFll 


87.7 


77.9 


246 ± 34.6 






BDFll 


342.8 


8.8 


865 ± 101 




387 ± 179 


BDFll 


42.5 


-11.8 


152 ± 30.7 






BDF12 


-241.0 


-256.7 


1620 ± 105 


-97.4 ± 368 




BDF12 


260.2 


300.5 


1191 ±129 


-1080 ± 905 




BDF12 


-137.4 


-133.4 


278 ± 40.5 


126 ± 36.7 


165 ± 273 


BDF12 


170.9 


66.7 


211 ±38.7 


98.8 ± 33.9 


-151 ± 259 


BDF13 


181.4 


-49.5 


721 ±51.2 


531 ± 50.8 


935 ± 271 


BDF13 


-154.0 


299.7 


1145 ±98.8 


126 ± 427 




BDF13 


225.1 


-99.9 


317 ±56.4 


11.4 ± 106 





Entries are left blank for fields not observed at 8 GHz. For the 30 GHz BIMA observations, entries are blank for those point 
sources which lie outside the primary beam and are not detected at > 3a significance. 
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TABLE 3: Point Sources and Fluxes Identified from VLA Observations(2002) 



Field 


A R.A. (") 


A DEC. (") 


4.8 GHz (/iJy) 8 GHz (^ijy) 30 GHz (jttJy) 


BDF14 


125.1 


-166.8 


1 1 86 + 35 3 


1 235 + 63 1 


BDF14 


-56.3 


30.3 


oofi + 24 1 


nnK + 22 2 — 1 95 + 1 43 


BDF14 


17.6 


-183.4 


578 + 96 5 


Qcy J. 4.Q q 1 1 n + 944 


BDF15 


-6.9 


-324.4 


q^qn + 78 fi 




BDF15 


-129.0 


218.7 


CQ7 -L 42 1 

yjKj 1 —1— ^£j.±. 




BDF15 


77.4 


-282.6 


707 ± 51.4 




BDF16 


-157.8 


116.7 


1091 ± 38.3 


212 ± 274 


BDF16 


-266.2 


198.5 


2553 ± 74.8 




BDF16 


-158.4 


33.6 


576 ± 33.1 


215 ± 222 


BDF17 


39.9 


169.3 


431 ± 26.2 


774 ± 44.4 694 ± 228 


BDF18 


-213.4 


381.6 


6902 ± 145 




BDF18 


-106.1 


-0.7 


465 ± 25.7 


145 ± 168 


BDF18 


285.4 


2.8 


988 ± 49.9 




BDF18 


137.7 


93.9 


515 ± 29.7 


103 ± 222 


BDF18 


-9.5 


290.3 


1636 ± 50.5 




BDF18 


112.9 


165.9 


266 ± 32.8 


381 ± 276 


BDF18 


- 32.2 


105.4 


333 ± 25.5 


142 ±170 


BDF19 


31.9 


169.4 


289 ± 32.8 


-95 ± 237 


BDF19 


-169.8 


- 17.3 


363 ± 32.6 


-16 ±234 


BDF19 


-128.7 


- 76.5 


403 ± 30.4 


167 ± 210 


BDF20 


245.8 


15.9 


3809 ± 49.4 


1240 ± 400 


BDF20 


-131.5 


394.5 


4957 ± 136 




BDF20 


89.7 


360.3 


3272 ± 96.6 




BDF20 


-301.9 


149.2 


2851 ± 76.0 




BDF21 


277.3 


-180.0 


2103 ± 57.5 




BDF21 


229.4 


74.1 


1187 ±34.8 




BDF21 


-315.9 


-224.9 


838 ±82.2 




BDF21 


115.7 


268.9 


470 ± 44.2 





Entries are left blank for fields not observed at 8 GHz. For the 30 GHz BIMA observations, entries are blank for those point 
sources which lie outside the primary beam and are not detected at > 3cr significance. 
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3. EXCESS POWER ESTIMATE 

Following the method of D2002, the excess power in 
the BIMA data is computed assuming that the angular 
power spectrum can be described by either one or two 
flat band powers. Point sources identified with the VLA 
are tabulated in H2.2l and are removed from the data us- 
ing the constraint matrix technique described in D2002. 
Confidence intervals for the band powers are determined 
using the integrated likelihood also described in D2002. 
In Tabled we show the most likely AT^ and approximate 
uncertainty computed from visibilities in the 0.63 — 1.1 kA 
range for each of the 18 fields in the survey. Results for 
the complete data set are included in TableOfor the cases 
when the power spectrum is modeled by two bandpowers 
corresponding to u-v ranges 0.63 — 1.1 kA and 1.1 — 1.7 kA, 
and for a single bin covering the range 0.63 — 1.7 kA. 

In the combined analysis of all 18 fields, we allow for 
AT^ < 0.0 in determining the most likely estimate of 
measured power and confidence intervals. Excess power 
corresponding to AT^ > 0, was observed with 94.7% 
confidence in the 0.63 — 1.1 kA bin. In the u-v range 1.1 — 
1.7 kA, the level of observed power was consistent with 
zero and we found an upper limit AT2 < 880/iK^ at 95% 
confidence. When an analysis is performed combining all 
data into a single bin, we find an estimate of excess power 
excluding zero, AT^ > 92.3% confidence. 

Window functions for each of these bands are produced 
from the noise weighted sum of the window functions for 
the individual visibilities. Averaged over all 18 fields, 
the 0.63 — 1.1 kA band has an average value of £eff = 
5237 with FWHM £ = 2870. The window function for 
visibilities in the u-v range 1.1 — 1.7 kA has an average 
value £eff = 8748 with FWHM £ = 4150. For the single 
band model covering the u-v range 0.63 — 1.7 kA, the 
window function has an average value £ef / = 6864 with 
FWHM £ = 6800. 

TABLE 4: Power Estimates for 0.63 - 1.1 kA 



Field 


Most Likely a 


BDF4 


0.0 


600 


HDF 


0.0 


270 


BDF6 


380 


455 


BDF7 


300 


1050 


BDF8 


0.0 


390 


BDF9 


0.0 


590 


BDFIO 


0.0 


360 


BDFll 


50 


920 


BDF12 


1590 


1165 


BDF13 


1480 


1260 


BDF14 


690 


920 


BDF15 


0.0 


660 


BDF16 


300 


880 


BDF17 


1390 


1400 


BDF18 


200 


990 


BDF19 


0.0 


900 


BDF20 


290 


910 


BDF21 


0.0 


310 



The likelihood distribution near the maximum (ATj = 
AT^s) for the data defined by bin B is well described by 
a offset log- normal function (Bond, Jaffe, & Knox, 2000), 

ln£(AT2) = ln£(AT^) - 1 ^ e^z. , 

B ^ 

(1) 

where the offset log-normal parameters Z are defined as 
ZB=\n{ATl+XB). (2) 

The likelihood functions are fit with this model to deter- 
mine values for the curvature at peak ub, which rep- 
resents the uncertainty in the measurement, and log- 
normal offset xb the likelihood functions reported in 
Table El 

3.1. Point Source Removal 

As discussed in t|2.2l we have adopted a detection 
threshold of 6cr for identifying point sources in the 
VLA data. We measured the effect of eight different 
point source detection thresholds on the measured ex- 
cess power from the combined analysis of the 18 fields. 
We choose point source detection limits in terms of SNR 
rather than flux to account for attenuation by the VLA 
primary beam. The noise is assumed to have an RMS of 
25/iJy/beam in all VLA fields. The results are listed in 
Table El 

In the case for which no point sources are removed, the 
most likely value of ATI = 430 ^lK^ appears significantly 
elevated due to contamination by radio point sources. 
The measured excess power drops from 430 /iK^ in the 
case of no point source constraints to a broad minimum 
of ^ 200 — 250 /^K^ for removal of sources detected with 
less than 8(7 significance. As the detection threshold de- 
creases, the removal of point sources will begin to remove 
a significant number of degrees of freedom from the anal- 
ysis, resulting in increased uncertainty in the measure- 
ment of power. For example, a detection threshold for 
point sources of Aa removes three times as many sources 
from the data as a detection threshold of 8 — 12a. The 
uncertainty in the measurement begins to significantly 
increase with a point source detection threshold less than 
6(7 while the subtraction of the additional sources have 
no effect on the observed excess power. Therefore, we 
adopt 6cr as the point source detection threshold. 

While the contribution to anisotropy from point 
sources is expected to scale as the estimate of ex- 
cess power is consistent with zero on finer angular scales 
as shown in Table El It should also be noted that no ad- 
ditional point sources were found in the 8 GHz VLA ob- 
servations. These observations reinforce the conclusion 
that the contribution of point sources to the observed 
excess power is well constrained in the power spectrum 
analysis. 

3.2. Systematic Tests 

We performed tests for systematic errors in all fields 
identified as having significant excess power. This anal- 
ysis was limited to visibility data in the u-v range 
0.63 — 1.1 kA described by AT^ where the most significant 
detection of excess power occurs. The modeled power in 
the second bin is fixed at Ar| = for all tests described 
in this section. We searched for systematic errors by re- 
peating the power spectrum analysis after splitting the 
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AT'^^K'-^) Likelihood 

u-v range Most Likely 68% Confidence as xb AT^ > 

0.63- 1.1 kA 220 100 - 360 130 625 94.7% 

l.l-1.7kA -40 < 420 395 3040 

0.63-1.7kA 170 70 - 290 110 490 92.3% 



TABLE 6: Effect of Point Source Model on IST^ 
VLA Detection Limit (SNR) Number of Sources Most likely cr 



none 





430 


155 


> 40 


7 


350 


150 


> 20 


27 


290 


140 


> 12 


45 


260 


130 


> 8 


58 


240 


135 


> 6 


62 


220 


130 


> 5 


98 


210 


135 


> 4 


168 


250 


150 



data into a number of subsets. The first set of tests 
was designed to search for systematic contamination that 
changes with time. Relative to a source on the celestial 
sphere, terrestrial sources will appear to move rapidly 
over the course of a single observation, while the sun 
or moon will vary in position by many degrees over the 
course of a typical month long observation. We searched 
for such signals by analyzing the data after dividing it 
into subsets corresponding to the first, second, and third 
sections of both observation tracks and the period of ob- 
servation. For the field BDF6, the only field observed 
in multiple years that was found to have a significant 
level of excess power, we also compared the results of 
observations taken in 1998 and 2000. Instrumental ef- 
fects that manifest themselves as spurious signals on a 
given telescope or baseline are also a potential source of 
systematic error. We searched for such effects by break- 
ing the data into subsets of four and five telescopes and 
looking for antenna based systematic errors. For a test 
of baseline based systematic errors, we created east- west 
and north-south baseline subsets. The different data and 
instrument subsets used in the search for systematic er- 
rors are listed in Table [7| The results for the application 
of these tests applied to fields BDF6, BDF12, BDF13, 
BDF14, and BDF17 are shown in FigurcUl 



TABLE 7: Cuts Used in Systematic Tests 



Test Number 


Data Subset 


1 


First Four Hours UT 


2 


Middle Three Hours UT 


3 


Final Four Hours UT 


4 


First Third of Observation Dates 


5 


Middle Third of Observation Dates 


6 


Final Third of Observation Dates 


7 


Baselines from Subarray of four Telescopes 


8 


Baselines from Subarray of five Telescopes 


9 


East- West Baselines 


10 


North-South Baselines 


11 


First Year of Observation (BDF6 only) 


12 


Second Year of Observation (BDF6 only) 
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Dato Subset ^ 
BDF13 



10 12 



Dato Subset ^ 




10 12 



Data Subset i 



Dote Subset ^ 



10 12 



Fig. 1. — : Band power estimates and confidence intervals for each systematic test in those fields found to have 
significant levels of excess power. The shaded region represents the 68% confidence intervals for the full data set of 
each field. 
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If a systematic error was associated with one of these 
subsets, then the level of excess power in that subset 
should increase. However, as can be seen in Figure^ we 
found approximately the same level of excess power in 
each subset. Of the fifty-two subsets, only the final four 
hours of UT in BDF6 (third data point in Figure ^) and 
the middle third of observation dates for BDF14 (fifth 
data point in Figure ^ are found to have estimates of 
excess power and 68% confidence intervals that lie above 
the 68% confidence limits of the full analysis. Therefore, 
we consider these two subsets to be the most likely to be 
systematically biased. To test the possible contribution 
from systematic errors in these subsets to the measured 
excess power, we repeated the analysis on the 18 com- 
bined fields after removing these two subsets. We found 
an estimate of AT^ = 180^12° mK^, not significantly dif- 
ferent than the reported value of AT^ ^ 220^120 for 
the entire survey. 

3.3. Correlated Signal Between Independent Fields 

It is possible that a hardware malfunction could create 
a systematic false correlation that is constant or changes 
slowly with sky position. To test for this, we combined 
the raw visibilities from observations of all fields taken 
in a single summer, as if all the data came from a sin- 
gle pointing. If the fields contain random noise or in- 
dependent sky signal, the power in the combined fields 
should decrease significantly when the visibilities are av- 
eraged. A false detection caused by correlations intro- 
duced in the hardware or other local effects might, de- 
pending on its stability, enhance the excess power in the 
anisotropy measurement when the independent observa- 
tions are combined. The results of the analysis of the 
combined data sets, listed in Table|Hl are consistent with 
instrumental noise at 68% confidence, as expected for 
non-correlated independent observations. 

To test for false correlations that vary slowly across 
the sky, we observed a set of fields offset by 2' in RA 
from fields that had been previously observed. The fields 
BDF6, BDF12, BDF13, BDF14, and BDF17 have the 
most significant levels of excess power of the individual 
fields in the survey. We have neighboring fields for all of 
these except for BDF13. The raw visibilities of the neigh- 
boring fields were combined and then analyzed for excess 
power; the results of this analysis are shown in Table |HI 
In all cases, the observed anisotropy power decreases as 
we would expect for the combination of visibilities from 
independent patches of sky. Therefore, we conclude that 
there does not appear to be any false correlation that 
is constant or slowly varying with sky position. In ad- 
dition, there does not appear to be any correlation be- 
tween excess power and proximity to the sun, none of the 
fields analyzed for systematic error were near the sun's 
location of 11 — 12 hr R.A. during the summer months. 
Overall, we find no evidence that our results are biased 
by systematic effects or astrophysical contamination. 

3.4. Constraints on erg 

As was described in the introduction, the most likely 
astrophysical source of the excess power observed in the 
survey is expected to be CMB anisotropy arising from 
the SZ effect in clusters of galaxies. Assuming that the 
observed excess power is entirely due to the SZ effect, we 



TABLE 8: Results of Combining Independent Observa- 
tions 
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use a publicly available archive of N-body simulations ^ 
(Schulz & White, 2003) to add the expected signal from 
the SZ effect to Monte Carlo noise realizations. These 
simulations predict the evolution of mass from z = 60 to 
present in a cosmological model described by VIm = 0.3, 
l^A = 0.7, n^K^ = 0.02, h = 0.7, n = 1 and CTg = 1.0. 
Only dark matter is included in these simulations, and 
the baryons contributing to the SZ signal are added as- 
suming that the gas closely traces the dark matter. In- 
dependent regions of the N-body simulations arc multi- 
plied by the 6.6' FWHM primary beam and transformed 
into the u-v plane with the same u-v sampling as the real 
BIMA data. Noise is added to each u-v point with a vari- 
ance determined from the observed visibilities. There- 
fore, each simulated observation has u-v coverage and 
noise characteristics identical to the real BIMA obser- 
vation of each field. The analysis of one hundred real- 
izations of the BIMA survey with unique instrumental 
noise and simulated SZ sky resulted in an excess power 
of AT2 = 216 ± 190 /iK^ at 68% confidence; this is re- 
markably close to the level of excess power found in the 
BIMA survey. 

Komatsu and Seljak (2002) demonstrate that the am- 
plitude of the SZ power spectrum has a strong depen- 
dence on (Tg, with AT^ oc cr| for erg near unity. We scale 

7 /2 

the simulated SZ images by tig to produce skies cor- 
responding to different values of erg. These images are 
transformed to the u — v plane, combined with the simu- 
lated instrumental noise, and used to compute the likeli- 
hood of each value of erg resulting in the observed excess 
power, which we approximate as /^T^ = 180 — 260/iK^. 
The relative likelihood that simulations with a given 
value of (jg will reproduce the observed excess power is 
determined from 10,000 realizations of the BIMA survey 
with tjg ranging from 0.0 to 1.5 The resulting relative 
likelihood shown in Figure |21 includes contributions to 
the uncertainty from both noise in the measurement and 
sample variance due the non Gaussian nature of the SZ 
signal and the small patch of sky surveyed. Assuming 
an additional 10% uncertainty in the simulations (Ko- 

^ Data available at |http: //pad. berkeley.edu/tSZ7] 
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matsu k Seljak, 2002, Goldstein et al. , 2003), we find 
0-8 = l-03to:29 at 68% confidence and as = l-03lo;96 
95% confidence. 



at 



Fig. 2. — : Relative likelihood of measuring the observed 
excess power as a function of as- The uncertainty at 
each point is determined by the Poisson statistics for the 
number of events in that bin. 
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FiG. 3.—: 

Typical synthesized beam created from the u-v coverage 
used in the analysis of image statistics. The dotted line 
shows the half-power radius of the primary beam. 

TABLE 9: Properties of Each Image Used in Analysis 
of Image Statistics 



4. IMAGE ANALYSIS 

Analysis of the complete BIMA survey results in a 
detection of excess power with nearly 95% confidence. 
However, with the power spectrum analysis, we are un- 
able to determine if the excess power is caused by sources 
with negative flux as would be expected from the SZ ef- 
fect in galaxy clusters at 30 GHz. The statistics of the 
pixel flux values in the images contain additional infor- 
mation that can in principle constrain possible sources 
of the excess power. 

Images are reconstructed directly from the visibil- 
ity data for each field in the BIMA survey using the 
DIFMAP software package (Pearson et al. 1994). Po- 
sitions determined from the VLA data are used to model 
and remove point sources (> 6a) in the u-v plane us- 
ing baselines with u-v radius > 1.5 kA (~ 50% of the 
visibilities). BIMA fields are imaged at a resolution of 
7.5" per pixel after applying a Gaussian taper with a 
half-power radius of 1.0 kA to the visibility data to max- 
imize brightness sensitivity. These mapping parameters 
are used when presenting images of bright SZ clusters 
observed with the BIMA array (i.e. Laroque et al. 2003) 
and produce a synthesized beam which is well matched 
to cluster scales. A map of a typical synthesized beam 
can be found in Figure The response of the first neg- 
ative sidelobe is a factor of four lower than response at 
the center of the beam. 

The typical RMS noise is 110/iJy/beam for an image 
with 90" FWHM synthesized beam. The RMS for each 
image listed in Table El are computed directly from the 
noise properties of the images. The RMS temperature 
estimates correspond to the synthesized beamsize com- 
puted for the u-v coverage of each experiment. 

To provide a visual representation of the data used in 
the analysis, images of the BIMA fields are reproduced 
in Figure^ The dashed circle in each image represents 
the radius at which the primary beam attenuates the sky 
signal by a factor of two relative to the pointing center. 
Regions lying far outside the primary beam can be used 
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to compute the level of instrumental noise in the map. 
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Fig. 4. — : BIMA Images after applying a Gaussian taper with a half-power radius of l.OkA to the visibility data. 
Point sources identified in the VLA images have been modeled from the long baselines {u-v > 1.5kA) and removed 
directly from the visibility data. Dashed lines correspond to the half power radius of the primary beam. 
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Fig. 4.—: Contd. 
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4.1. Image Simulations 

Decrements caused by galaxy clusters would be ex- 
pected to produce an excess of high SNR negative pixels. 
We attempt to detect this unique signature using higher 
order statistics of the image flux distribution as described 
in Rubino-Martin & Sunyaev, 2003. We compute asym- 
metry, skewness, and extrema in the observed data set. 
These results are compared to both simulations of in- 
strumental noise and to simulations of SZ galaxy clus- 
ters added to instrumental noise in order to determine 
the significance of the observed distribution of flux. 

The Fourier transform of the visibilities measured by 
the interferometer into the image plane introduces corre- 
lations between pixels which complicate the noise prop- 
erties. In order to better understand this, we generated 
Monte Carlo simulations of the visibility data and trans- 
formed them into images. We created two sets of sim- 
ulations, each containing 100 realizations of the full set 
of 18 images. These simulations are used to quantify the 
significance of applying the statistical tests to the BIMA 
images. The first set of simulations uses random complex 
visibilities with variances consistent with the observed 
noise for each image. At each point in u-v space in the 
observed data, a simulated visibility is created with a 
real and imaginary component from a random sampling 
of a Gaussian distribution with variance determined from 
the weight of that visibility. Therefore, each simulated 
observation has u-v coverage and noise characteristics 
identical to the real BIMA observation of each field. 

In the second case, we add the expected contribution 
from the SZ sky to the MC noise realizations using the 
same simulated SZ sky images described in Section 
with (Tg fixed to 1.0. Independent regions of the N-body 
simulations are attenuated with the 6.6' FWHM primary 
beam and transformed into the u-v plane using the exact 
same u-v sampling as that in the real BIMA data. Noise 
is then added to the visibility data in the same way as 
for the noise-only simulations and the data is then trans- 
formed into the image plane for analysis. 

The image statistics of the observed images are com- 
pared with those generated from the Monte Carlo simu- 
lations. For each simulation, an image is generated from 
the raw visibilities tapered in the u-v plane with a 1.0 kA 
FWHM Gaussian. Statistics are generated directly from 
the images of the 100 simulations. Using a large number 
of Monte Carlo simulations, one can produce a probabil- 
ity distribution for each observable statistic. The most 
likely value and integrated confidence intervals can be 
taken from this distribution and compared to the results 
of the real observations. Due to the modest number of 
100 simulations of the survey, the shape of the pixel dis- 
tribution is strongly dependent on bin size and bin spac- 
ing, leaving the most likely value as a poorly determined 
quantity. Instead, the median value in the 100 simula- 
tions is taken to be the most likely estimate. The terms 
defining the 16th and 84th percentile in the distribution 
are taken to be the lower and upper bounds to the 68% 
confidence intervals, respectively. The 4 nearest neigh- 
bors to the median, lower bound, and upper bound in 
the distribution are averaged to reduce the noise in the 
estimates. 

We use a histogram of pixels binned in intervals of 
significance, S = Xi/a, where Xi is the value at pixel 



i and a is the estimated image RMS (approximately 
110/zJy/beam). This is equivalent to weighting the pix- 
els by the sensitivity of the observation. We can compare 
the pixel distributions for observed data and simulations 
by comparing the histogram of the images as shown in 
Figure 13 
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Simulations of Noise in BIMA Survey 




Simulations of Clusters in BIMA Survey 




Fig. 5.—: 

The figure on the left represents the histogram for simulations including only instrumental noise. The figure on the 
right represents the simulations with clusters added. In both cases, the solid line represents observed data, the dashed 
line represents the most likely value from the simulations, and the dotted lines represent the 68% confidence interval 
of the simulations. The simulations including clusters are a better fit to the observed data not only in the region of 
negative pixels, where the SZE is expected to directly contribute to the signal, but also in the region of positive pixels 
where the effects of the negative sidelobes contribute to the signal. 
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4.2. Asymmetry of Pixel Distribution 

For a given image, we can characterize the pixel flux 
distribution by selecting a flux interval AD, and com- 
puting a histogram of the number of pixels with a flux 
between D - AD/2 and D + AD/2). The asymmetry 
of this histogram can be estimated directly as the differ- 
ence in area between the positive and negative regions. 
A comparison of this statistic for the BIMA images and 
simulations is shown in Figure 

Simulations of BIMA Survey 




-200 r I I - 

-6 -5 -4 -3 

Fig. 6.—: 

This figure describes the asymmetry of the pixel distri- 
bution as a function of SNR. The dark shaded region 
and dotted line represent the 68% confidence intervals 
and median asymmetry of the simulations that include 
only instrumental noise. The region of angled parallel 
lines and the dashed line represent the 68% confidence 
intervals and median asymmetry of simulations that also 
include the SZE. The black solid line represents the ob- 
served data. 

It should be noted that the cluster simulations produce 
an excess of positive pixels compared to simulations of in- 
strumental noise only. There are no sources with positive 
flux included in the cluster simulations, however pixels 
with positive flux are expected due to the sidelobes of 
the synthesized beam. A Gaussian taper was applied to 
the u-v data to minimize the sidelobes and this effect. 

It is clear from Figure El that the observed data is con- 
sistent with the asymmetry in the simulations that in- 
clude galaxy clusters and inconsistent with simulations 
that include only noise. For example, the excess of pix- 
els with S" < — 3 in the observed data exceeds that in 
the simulations of instrumental noise at 86% confidence. 
The number of pixels with S < —4 exceeds that found in 
the simulations of only instrumental noise at 88% confi- 
dence. Therefore, the application of this statistic results 
in a detection of signal with the morphology of the SZ 
effect with a significance comparable to the significance 
of the detection of excess power reported in Section|31 In 
fact, the observed asymmetry in the BIMA data exceeds 
the mean value from the noise & SZ simulations. 

As a simple exercise, a similar analysis is used to rule 
out positive flux (such as point sources) as the cause 
of excess power with slightly less confidence; we simply 
repeat the analysis described, but with the sign of the 
flux from clusters reversed. The excess of pixels with 
5* < — 3 in the observed data exceeds that in the simula- 
tions using positive flux at 74% confidence. The excess 



of pixels with S < —4 exceeds that in simulations of 
instrumental noise at 76% confidence. A more rigorous 
analysis for point sources would include a more realistic 
distribution of point sources and take into account the 
details of their removal using the full, equally weighted, 
set of visibilities. This expanded analysis is considered 
unnecessary given the results of the VLA observations, 
constraints on point sources, and lack of power at finer 
angular scales. Nonetheless, this method would be useful 
in cases where rejection of point sources was less certain. 
We again conclude that it is very unlikely that point 
sources are responsible for the asymmetry observed in 
the BIMA image pixel fluxes. 

4.3. Skewness of Pixel Distribution 

A measure of skewness conveys information about the 
sign of the features producing the deviation from Gaus- 
sianity. This quantity can be estimated using the third 
moment of the data: 

The skewness of the data derived from the BIMA images 
is compared to the skewness determined from the Monte 
Carlo simulations in Figured 

Simulations of DIMA Survey 
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Fig. 7.—: 

The cumulative count of simulations as a function of 
skewness. The solid line represents the model assuming 
only instrumental noise and the dashed line represents 
the model including the simulated SZ sky from the n- 
body simulations. The vertical line represents the skew- 
ness measured in the BIMA data. 

Analysis of the simulations of instrumental noise alone 
results in a skewness of 0.01 ± 0.08 at 68% confidence, 
while the simulations of clusters plus instrumental noise 
give a mean skewness of — 0.096lg ijg. The skewness in 
the BIMA observations was found to be F = —0.066, 
a value that is inconsistent with the noise only simula- 
tions at 64% confidence. Repeating the simple analysis 
described at the end of ii4.2l the skewness is inconsistent 
with positive flux sources as the cause of the observed 
excess power at 84% confidence. 

4.4. Outliers in the Distribution 

In addition to the two tests described above, we per- 
formed simulations to characterize the significance of out- 
liers in pixel flux distributions of the individual fields. 
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We identified the most negative pixels in each of the 18 
BIMA fields and compared these with the results of the 
simulations. Seven BIMA fields were observed to have 
decrements with S < —3a, more than occurred in 93% 
of the survey simulations with noise only. Two BIMA 
fields were observed to have decrements with S < —Aa, 
more than occurred 99% of the survey simulations in- 
cluding only noise. None of the observed BIMA fields 
had a pixel with positive flux S > 4tT. 

4.5. Interpretation of Image Statistics 

In each of the statistical tests described above, the 
BIMA images were found to be inconsistent with those 
produced by instrumental noise alone at approximately 
the l — 2a level, comparable to the significance of the de- 
tection of excess power. Although the analysis does not 
conclusively determine the source of the excess power, 
the results are consistent with the signal expected from 
SZ galaxy clusters. The results are also inconsistent with 
radio point sources being the source of the observed ex- 
cess power. Despite the lack of a definitive conclusion 
from these tests, the results are encouraging. In order 
to investigate what observations would be required to 
make a definitive measurement, we performed simula- 
tions of the 18 field survey with four times the obser- 
vation time, or equivalently four times the correlation 
bandwidth. These simulations resulted in a detection of 
excess power and skewness at the 4cr level in more than 
50% of the simulations. Such a detection would provide 
convincing evidence for SZ clusters as the source of excess 
power. 

5. OPTICAL OBSERVATIONS OF BIMA FIELDS 

To date, all known galaxy clusters have been discov- 
ered through optical or X-ray observations. A negative 
correlation of the observed fine scale CMB anisotropy 
with X-ray or optical emission would be a smoking gun 
for the discovery of a cluster through the SZ effect. To 
search for this correlation, we observed ten of the BIMA 
fields using ground-based optical telescopes. We selected 
the five fields with significant levels of excess power as- 
suming that they are the most likely candidates for iden- 
tifying galaxy clusters. In addition to those five fields, 
we observed five BIMA fields that lie at convenient RA 
during the nights we were awarded time. Imaging was 
performed using I and R filters on the LRIS instrument 
(Oke et al. 1995) on the 10 m Keck^ I telescope. Imag- 
ing in z' was done with the MOSAIC instrument (Wolfe 
et al. 1998) on the 4m Kitt Peak National Observatory 
telescope^. 

We performed an analysis of the optical images us- 
ing a method similar to that used in the Red Cluster 
Survey (RCS) survey (Gladders & Yee, 2000). We first 
categorized objects into two redshift bins determined by 
color R — I and R — z' . The first bin contains only high 

® Keck Observatory is operated as a scientific partnership among 
the California Institute of Technology, the University of California 
and the National Aeronautics and Space Administration. The Ob- 
servatory was made possible by the generous financial support of 
the W.M. Keck Foundation. 

^ KPNO is a Division of the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc. (AURA) under cooperative agreement 
with the National Science Foundation. 



redshift galaxy candidates. These candidates have color 
i? — / > 1.0 and i? — z' > 1.0, implying a red sequence 
redshift 0.5 < z < 1.0. The second bin contains low 
redshift galaxy candidates. These candidates have color 
0.0 < i? - / < 1.0 and 0.0 < i? - z' < 1.0 implying a red- 
shift 2 < 0.5. All objects not satisfying these criteria are 
considered foreground contamination and are discarded. 
An exponential kernel is used to smooth the maps of de- 
tected objects in each redshift bin with scale radius of 
20" to create a surface density map. We then create a 
product map by simply multiplying each surface density 
map with the corresponding BIMA map described in 21 
The product maps are searched for peaks and asymmetry 
in the distribution. Statistics are quantified with Monte 
Carlo simulations of random BIMA fields. The result- 
ing analysis (described in full detail in Dawson, 2004) 
showed no significant correlation between overdensities 
of galaxies in the optical maps and decrements of flux in 
the BIMA 28.5 GHz maps. 

While the follow-up optical observations do not con- 
firm the hypothesis that the observed anisotropy is 
caused by the SZ effect in galaxy clusters, the method 
of combining optical and SZ observations should prove 
to be a powerful technique for identifying galaxy clusters 
in future SZ surveys. In this first attempt at doing so, 
we have not modeled the expected optical cluster signa- 
ture sufficiently to say what constraints this null result 
places on the role of SZ clusters in producing the ob- 
served excess power. There also remains the possibility 
that anisotropy from galaxy clusters at redshifts beyond 
the sensitivity of the optical data, z > 1, could be con- 
tributing significantly to the observed signal in the BIMA 
survey. 

6. CONCLUSION 

In this paper, we report the final results from our 
search for arcminute scale CMB anisotropy using the 
BIMA array. Modeling the observed power spectrum 
with a single flat band power with average multipole 
of = 6864, we find AT^ = nQtf^^lK^ at 68% 
confidence and a detection of AT^ > at 92.3% con- 
fidence. Dividing the data into two bins corresponding 
to different spatial resolutions in the power spectrum, 
we find ATf = 22{)t\^^JK^ at 68% confidence for CMB 
flat band power described by an average multipole of 
l^ff = 5237 and Ar| < MO ^iK^ at 95% confidence 
for ieff = 8748. We have used VLA observations and 
various cuts to test for contamination from radio point 
sources and systematic effects and conclude that it is un- 
likely that these sources are responsible for the observed 
signal. If we assume that the measured excess power is 
due to a background of distant SZ clusters, we can com- 
pare its value with that from simulations of large scale 
structure to place a constraint on the normalization of 
matter fluctuations, erg = 1.03lQ'2g ^8% confldence. 

In order to try to determine the source of the ob- 
served anisotropy power, we have performed an analysis 
of the BIMA image statistics. We compared the skew- 
ness, asymmetry, and outliers of the measured pixel flux 
distribution with simulations including noise only and 
noise plus SZ clusters. A statistical analysis of the BIMA 
survey images found that they were consistent with sim- 
ulations including a background of SZ clusters, and in- 
consistent with simulations of instrumental noise alone or 
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noise plus radio point sources at 1— 2(t. Additional Monte 
Carlo simulations indicate that with approximately four 
times the time dedicated to the survey, or equivalently 
four times the correlated bandwidth, the BIMA instru- 
ment would achieve the sensitivity to test the hypothesis 
of SZ clusters as the source of the observed excess power 
at greater than 99% confidence. Therefore, future ded- 
icated interferometers, such as the Sunyaev-Zel'dovich 
Array-^'^, should be able to effectively use image statis- 
tics to determine if any observed anisotropy is due to 
SZ clusters. Finally, we performed a preliminary search 
for a correlation between red galaxy density and CMB 
temperature fluctuations. We are currently unable to 
quantify the significance of the null result, however, we 
expect that X-ray and optical follow-up will be essential 
tools for the interpretation of future SZ surveys. 

We thank the entire staff of the BIMA observatory 
http://astro.uchicago.edu/sza/ 
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